Dust absorption along the line of sight for high-redshift objects 



M. Trent i 

and 
M. Stiavelli 

Space Telescope Science Institute, 3700 San Martin Drive Baltimore MD 21218 USA 
trentiOstsci . edu; mstiavel@stsci . edu 

ABSTRACT 

We estimate the optical depth distribution of dust present in absorption sys- 
tems along the line of sight of high redshift galaxies and the resulting reddening. 
We characterize the probability distribution of the transmission to a given red- 
shift and the shape of the effective mean extinction law by means of analytical 
estimates and Monte Carlo simulations. We present our results in a format use- 
ful for applications to present samples of high redshift galaxies and discuss the 
implications for observations with the James Webb Space Telescope. Our most 
realistic model takes into account the metallicity evolution of Damped Lyman a 
absorbers and predicts that the effects of dust absorption are modest: at redshift 
z > 5 the transmission is above 0.8 at an emitted wavelength A e = 0.14/im with 
probability 90%. Therefore dust obscuration along the line of sight will affect 
only marginally observations at very high redshift. 

Subject headings: dust, extinction - - galaxies: high-redshift - - intergalactic 
medium — galaxies: ISM 



1. Introduction 

Dust along the line of sights affects the observations of distant astronomical objects and 
its effects on both extinction and reddening need to be taken into account. In the case of 
high redshift quasars this problem was initially addressed by Ostriker & Heisler (1984) and 
it has been the subject of several investigations thereafter. Fall & Pei (1989, 1993) developed 
a comprehensive theoretical framework to quantify the effects of dust along the line of sight 
and characterize sample selection effects induced by Damped Lyman a (hereafter DLA) 
absorbers. The reddening measurements by Pei et al. (1991), based on a sample of 13 DLA 
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systems with average absorption redshift (z a b s ) ~ 2.6, evidenced a reddening of background 
quasars with DLA spectral fingerprints with respect to a control sample without intervening 
absorbers. The reddening was measured from the shift of the average slope of the quasar 
spectral energy distribution between the two samples, which was found to be (Aa) = 0.5 
(significant at above 99% of confidence level). 

Some recent determinations based on larger samples of quasars (Murphy & Liske 2004; 
Ellison et al. 2005) do not confirm the earlier conclusions by Pei et al. (1991) and the new 
limits set on the dust obscuration along the line of sight are down by one order of magnitude 
with respect to the earlier estimates. Murphy & Liske (2004) find that Aa < 0.2 at 3 o for a 
larger sample of absorbers at (z a b s ) ~ 3. This implies a reddening E(B — V) < 0.02mag (at 
3a) . Unfortunately, the determination of the dust absorption at high redshift is intrinsically 
an indirect measurement and there is the possibility that the real dust content is higher than 
these lower estimates due to observational biases such as sample selection effects. Indeed, for 
a subsample of Damped Lyman a absorbers, selected by the presence of the Call absorption 
line, Wild & Hewett (2005) and Wild et al. (2006) find a significant evidence of reddening 
at moderate redshift {(E(B — V)) > 0.1 at z a b s ~ 1); similarly York et al. (2006) measure 
(E(B — V)) up to 0.085 for Mgll absorbers at z a b s ~ 1.4. Thus, any modeling of dust effects 
on high redshift objects will need to take into account possible selection effects. 

Surprisingly, very little consideration is generally given to the effects of dust obscuration 
along the line of sight for high-redshift, non- active galaxies. Clearly, if the effects of dust 
absorption have been detected for quasars, it is likely that they will affect every other object 
at similar distances. Nonetheless, it is common practice in observations of high redshift 
galaxies to consider the obscuration of dust as a screen localized at the emitter location 
(e.g., see Papovich et al. 2001), adopting a description of the dust properties like the one 
used for local starburst galaxies (Calzetti et al. 1994). 

The primary goal of this paper is therefore to estimate the fraction of essentially unob- 
scured lines of sight for very high redshift objects (z < 20) and their average transmission. 
The estimation of this quantity is much more robust with respect to uncertainties and biases 
in the observed distribution of absorbers than the measure of the amount of dust in DLA sys- 
tems (whose determination is outside the scope of this work). In fact, while a small number 
of optically thick absorbers, missed in magnitude limited surveys, could in principle contain 
the majority of dust in the universe, their effect on the average transmission along a random 
line of sight would be limited by their covering factor. The CORALS radio selected survey 
(Ellison et al. 2001) probes 66 lines of sight with complete optical follow-up detection. The 
probability of finding an optically thick absorber along a random line of sight is therefore 
below 4.9% at 99% of confidence level (and below 2.9% at 95% of confidence level). 
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In Section 2 we present our model for the absorbers, that is calibrated in Section 3 
up to redshift z ~ 5 on the measurements from recent Damped Lyman a (hereafter DLA) 
surveys (Ellison et al. 2001; Prochaska et al. 2005). For extrapolation to higher redshift we 
assume an unevolving comoving density of Damped Lyman a systems and a dust to gas ratio 
decreasing exponentially with redshift (Section 3.1). We also consider a wider set of input 
parameters to investigate different parameters extrapolation recipes and to study a higher 
dust content that could be missed due to selection effects. In Section 4, we describe our 
Monte Carlo code, that allows to characterize the full probability distribution of absorption 
which is presented in Section 5. Section 6 sums up. 



Given a source at redshift z e , we are interested in modeling the absorption due to dust 
residing in DLA systems (i.e. absorption systems with neutral hydrogen column density, N d , 
above 2 • 10 20 cm~ 2 ). Our approach is inspired by the models by M0ller & Jackobsen (1990) 
and by Madau (1995) (see also Fall & Pei 1993), in which the absorber distribution is treated 
as an input parameter that we calibrate to observations in the next section. 

We assume a discrete probability distribution of absorbers along the line of sight up to 
redshift z e with a separable probability distribution in column density ( N d ) and redshift (z) . 
This probability distribution is further assumed to be Poissonian per unit redshift, so that 



with E[dip(z)/dz] = n(z) (in this paper the symbol E p [f] means the expectation value of / 
under the probability measure p, and we drop the subscript p if it is clear what probability 
measure we are referring to). 

We consider the following forms for the column density distribution (f)(N d ): (i) a gamma 
function, the best fitting functional form observationally identified by Prochaska et al. (2005): 



with Nd in [N m i n ; +oo] and (ii) a power law (computationally friendlier in our Monte Carlo 
approach due to the simple analytical primitive) in a range [N min , N max \: 



2. Absorber modeling 



p(N d ,z)=(j>(N d )^(z), 



(1) 




(2) 



<f>(N d ) oc (N d ) 



(3) 



The normalization for the function 4>(N d ) is chosen so as to have J dN d Q(N d ) = 1. 
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The average number of absorbers per unit redshift n(z) is assumed to vary as: 

n(z) = AdX/dz, (4) 

where 

dx - H ° (1 1 z? - (i + *) 2 _ (5) 

Here we adopt a WMAP concordance cosmology with Q A = 0.7, Qm — 0.3 and H = 
70 km/s/Mpc (Spergel et al. 2006). 

The absorption cross section of dust located at z a is assumed to be a(X a ) in its rest 
frame. For an observed wavelength A G this can be written as: <r(A /(l + z a )). In the B-band 
in the absorber rest frame, A a = 0.44/im, we have that the optical depth of the dust can be 
written as: 

t b = N d a(0M»m) = k - , (6) 

10 Zi cm A 

where A; is a dimensionless dust to gas ratio parameter [k is of order unity for the Milky 
Way). 

At different wavelengths the absorption cross section can be expressed via an extinction 
curve £(A): 

£(A) = a(X)/a(X B ). (7) 
Thus, for a single cloud at redshift z a we can write the optical depth r(A G ) as: 

r(Ao) = io^^ (Ao/(1 + ^ )) - (8) 

Given these assumptions, we can compute analytically - down to the numerical eval- 
uation of a single integral over redshift - the average value for the transmission coefficient 
q(\ ,z e ) = exp {— r(A D , z e )} to a source at redshift z e . The computation is straightforward 
(see Appendix A in M0ller & Jackobsen 1990) and yields: 



E[q(X , z e )} = exp j- jT dz n{z){\ - ^[g(A G , z)\) 



(9) 



with 

q(X , z) = exp {-k{z)N d Z(X /(l + z))}. (10) 

We recall that E^[f] denotes the expectation value of / under the probability distribution 
0. 
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From the average value of the transmission coefficient we can define an effective optical 
depth: 



This quantity is a measure of the average departure from unit transmission and is more 
physically relevant than the mean optical depth for the purpose of characterizing the fraction 
of obscured lines of sight in the sky (see also Madau 1995). In fact the mean value of the 
optical depth is very sensitive to the detailed properties of the high tail in the r distribution. 
The concept can be easily illustrated with the following example. Consider 100 lines of sight, 
one of which with a really optically thick absorber (r = 10 4 ) while the others are optically 
transparent. In this case, the effective optical depth (Eq. 11) is 0.01 and correctly 

captures the fact that with probability 1% a line of sight is optically thick. The average 
optical depth is instead E[t] = 10, a very misleading value if applied to the estimation of 
the probability of having a line of sight free from absorption. 



We can take advantage of recent surveys (Ellison et al. 2001; Prochaska et al. 2005; 
Akerman et al. 2005; Rao el al. 2006), that have measured the observed redshift-column 
density distribution and metallicity for DLA systems, to critically examine the assumptions 
adopted in our model and to estimate its free parameters. 

In principle, the calibration of our model for the absorbers appears straightforward, as 
it relies on the observed statistics of Lyman a features in the spectra of several thousands 
of quasars over an extended redshift range. Unfortunately, for optically selected quasars 
samples, the observed distribution of DLA systems is in general a biased estimator of the 
intrinsic distribution as highly obscured lines of sight are preferentially missed for an optically 
selected sample. However, as we are interested in characterizing the average transmission to 
high z, and not in the measure of the comoving gas density in DLA systems, the dust bias at 
the high end of the column density distribution of neutral gas is of limited impact. In fact, 
if a small fraction e of optically thick absorbers is missed in a survey, this will only introduce 
a relative error of order e in the average transmission. The comoving gas density may on the 
other side be affected by an arbitrarily large error if these missed absorbers dominate the 
gas density budget. 

For the calibration we resort to two samples of DLA data. The first is the radio selected 
CORALS survey (Ellison et al. 2001) that has the advantage of being free from dust bias, 
but consists only of 66 quasars with detection of 19 intervening DLA systems. The statistical 




(11) 



3. 



Parameter estimation from DLA data 
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uncertainties in the values of the estimated parameters are rather large (reported as entries 
"E101" in Table 1). We therefore also consider the significantly larger, but optically selected, 
SDSS DR3 DLA dataset (Prochaska et al. 2005) that consists of 525 DLA systems identi- 
fied in the spectra of 4568 quasars with signal to noise ratio above 4. This dataset, while 
providing an excellent statistical accuracy (the observed neutral gas density in DLA systems 
is measured with relative error below 10%; see Prochaska et al. 2005) may be affected by 
systematic uncertainties for the high column density tail of the DLA systems distribution 
(see Trenti & Stiavelli 2006 for a characterization of the systematic errors in the SDSS DR3 
DLA dataset). The sets of parameters estimated using these data are reported as entries 
"Pr05" (SDSS) in Table 1. 

The absorbers column density distribution <&(N d ) has been characterized starting from 
the empirical distributions for the CORAL and SDSS surveys using a maximum likelihood 
estimator and plotting, in Fig. 1, the likelihood function. 

We estimate the values for the parameters in the gamma function description for (f)(N d ) 
using only the SDSS data, as the fit would have too many free parameters for the size of 
the CORALS survey. With our analysis we re-derive the same parameters identified by 
Prochaska et al. (2005) (entry "Pro05_T" in Table 1). Namely, we assume the standard 
DLA limit N min = 2 • 10 20 cm~ 2 , and we find a x = 1.8 and A" 7 = 3 • 10 22 cm~ 2 . For the 
single power law description we obtain a<i = 2.2 (SDSS) and a 2 = 2.1 (CORALS) adopting 
N m i n = 2 • 10 20 cm~ 2 (the standard DLA limit density) and N max = 10 22 cm~ 2 . This cut-off 
has been introduced with the goal of eliminating unphysical high N d tails in the distribution 
of gas column densities. The cut-off has been set to a value marginally higher than all the 
N d measurements in the SDSS and CORALS survey. We also explore different power law 
models with increasingly high value for the cut-off (models E101_a, E101_d and E101_e ) in 
order to assess the effects of a small additional number of optically missed absorbers with 
increasingly high column densities. While the SDSS data rule out a single power law at a 
confidence level above 3a for the observed column density distribution, this description may 
still be valid for the intrinsic distribution. Unfortunately, the CORALS data do not allow to 
significantly constrain the functional form of <fi{N d ) (see Ellison et al. 2001). The bias in an 
optically selected DLA survey maps an intrinsic power law distribution of column densities 
into an observed gamma function distribution (Fall & Pei 1993). 

In our approach, A is assumed to be independent of redshift, i.e. we are assuming 
a constant comoving number of absorbers. Due to the separability of Eq. 1, this implies 
a constant comoving density of neutral gas in DLA systems (f^/ A) ). The SDSS data 
(Prochaska et al. 2005) show evolution of by about a factor 2, mainly in the redshift 

range [2.2; 3]. However, at lower redshift, the measured gas density seems to be more in line 
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with the z > 3 values (see Fig. 22 in Prochaska et al. 2005 and Fig. 16 in Rao el al. 2006), so 
that this schematic modeling appears in reasonable agreement with the observations. The 
combined evolution seen by Rao el al. (2006) in the line and column density distributions 
at approximately constant Q^^ A ^ would affect the separability assumption of Eq. 1. This 
would complicate the numerical treatment of our model, but would not significantly change 
the value of the effective extinction to a given redshift, that, under the assumption of optically 
thin absorbers, depends in first approximation only on the comoving dust (i.e. neutral gas) 
density. The average transmission would depend more critically on the precise form of Eq. 1 if 
a significant population of optically thick absorbers is present, but this is an unlikely scenario 
given the complete optical follow-up detection in the radio selected CORALS survey Ellison 
et al. (2001). 

In order to estimate A, we have tested the covering factor of DLA systems. Starting 
from the published data we have identified for each quasar in the two samples the range of 
redshift in which the presence of DLA systems was detectable (from Table 1 in Prochaska 
et al. 2005 and from Table 3 in Ellison et al. 2001), computed the expected number of DLA 
systems in that interval and obtained A by evaluating the likelihood of getting the observed 
number distribution of DLA counts. 

The likelihood for A given the two datasets is reported in Fig 1: the maximum value 
for CORALS is at A = 0.0910, while the maximum for SDSS is at A = 0.0715. Assuming 
the SLOAN value for A, the CORALS data are marginally consistent: the higher CORALS 
result may be due to small number fluctuations at the la level (see the likelihood ratio in 
Fig. 1). An alternative possibility is that the discrepancy is due to an obscuration bias for 
SDSS. This is however highly unlikely, as the covering factor determination is dominated by 
DLA systems at the low end of the column density distribution, where the obscuration bias 
is negligible. The A value estimated from the SDSS data could even be an overestimate for 
the intrinsic A because of a Malmquist bias. That is, more absorbers with column density 
below the DLA limit have been scattered into the DLA sample than absorbers above the 
limit have been scattered out (O'Meara et al. 2006, Prochaska: private communication). 

For these parameters, we can compute the comoving density of neutral gas in DLA 
systems implied by our model (e.g. see Prochaska et al. 2005): 

CPc 

where mu is the mass of the hydrogen atom, \i = 1.3 is a correction factor for the composition 
of the gas, c the speed of light and p c the critical density of the universe. We have: n^ T LA> = 
0.81 • 10~ 3 for SDSS fitted to a gamma function, fi^ Lj4) = 0.84 • 10~ 3 for the SDSS fitted 
to a power law and = 1.17 • 10~ 3 for CORALS (with A = 0.0910). The agreement 
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with the published value from SDSS is excellent (see Table 9 in Prochaska et al. 2005: their 
unbinned measurement is (0.817 ±0.05) • 10~ 3 ), while our data cannot be compared directly 
with the analysis in Ellison et al. (2001), as they have used a different cosmology (0 M = 1 
and f^A = 0). 

The value of the dust to gas ratio k for DLA systems depends on their metallicity Z and 
as a first approximation we can consider a linear dependence of the dust to gas ratio on Z. 
DLA systems are generally characterized by a low metallicity and by a moderate evolution 
of their properties with the redshift (Wolfe et al. 2005). Their metallicity has been measured 
in several surveys and the average metallicity in optically and radio selected samples appears 
consistent (Akerman et al. 2005). Here we consider the compilation by Kulkarni et al. (2005) 
(see also Prochaska et al. 2003) and we approximate the reported measurements (Fig. 13 in 
Kulkarni et al. 2005) with a linear function for log (Z(z)), which provides a good agreement 
with the data in the redshift range 2 < z < 5: 

Z(z)/Z Q = 0.2 • lO" - 22 . (13) 

By considering a typical dust to gas ratio k = 0.8 for our galaxy (Z « Z Q , with ps 50% of 
the metals locked in dust grains), this translates into an observed dust to gas ratio: 

k(z) = 0.16-10~ - 22 . (14) 

To account for uncertainties in the measure of Z(z) and for a smaller fraction of metals 
depleted into dust, we introduce a free factor a K for the intrinsic dust to gas ratio k(z): 

k(z) = a K ■ k (z). (15) 

Realistic values for a K range in the interval [0; 1]. a K = 0.5 corresponds to 25% of the total 
metal amount in dust grains (e.g., see Pettini et al. 1997, Prochaska & Wolfe 2001), while 
a K — 1 implies a depletion factor like in the Milky Way (Pei & Fall 1995; Pei et al. 1999). 
Our main results are presented in terms of a K — 1, so that we effectively obtain upper limits 
on the obscuration along the line of sight. In the following Sections we will also discuss 
scenarios with a K < 1. 

The extinction curve £(A), in the absorber rest-frame, is assumed to be as measured 
in and parameterized for the Small Magellanic Cloud (Pei 1992). In fact, extinction curves 
similar to those measured for our Galaxy and for the Large Magellanic Cloud seem to be 
ruled out by the current observational data (Ellison et al. 2005; York et al. 2006); extinction 
curves for DLA systems at low redshift start to be directly measured and highlight a rather 
complex picture: Junkkarinen et al. (2004) have measured £(A) for a z = 0.524 DLA absorber 
finding some similarities with the Galactic extinction curve. For our purposes the use of a 
different extinction curve would not affect our results significantly, as the redshift averaging 
process over many absorbers tends to smooth out the specific features of the input curve. 
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3.1. Extrapolation of parameters up to z « 20 

As we are mainly interested in characterizing the expected absorption for future obser- 
vations at z > 6, the parameters that have been tuned to the properties of DLA systems at 
z < 5 have to be extrapolated into a redshift region with no observational constraints. 

Qualitatively, a monotonic metal (and dust) abundance appears plausible even before 
reionization when all hydrogen is neutral. Indeed, as the redshift increases the metallicity 
decreases, so one expects that the average local content of dust will progressively be reduced. 

Our derivation of the average transmission depends on the comoving dust distribution 
which we treat as the product of dust-to-gas ratio times neutral hydrogen distribution in 
discrete systems which we identify as DLA systems at z < 5 but that could simply be the 
sites of metal production at higher z. 

At a sufficiently large z the average comoving density of neutral gas will stay approxi- 
mately constant: the star formation rate should drop after z ~ 6 and this, combined with 
progressively less time available for star formation, means that the initial reserve of gas 
should remain almost undepleted. Eventually, in a hierarchical formation scenario, one can 
expect that the neutral gas will reside in more numerous Lyman a systems with smaller 
column densities, but this will influence only marginally the value of the expected absorp- 
tion along the line of sight at constant comoving gas density. In fact, given their very low 
dust-to-gas ratio, these primordial absorbers are expected to be optically thin, so that the 
average transmission will not depend in first approximation on the details of the distribution 
and will be proportional to the integrated comoving dust density. 

These arguments led us to choose to extrapolate our set of parameters with minimal 
assumptions: constant A, i.e. constant comoving VLhi density, and exponentially decreasing 
metallicity, as given by extension of our Eq. 14. These extrapolations are consistent with 
the asymptotic behavior of global models for the chemical evolution of DLA systems, like 
the ones developed by Pei & Fall (1995). 

The use of Eq. 13 implies a metallicity of 1.2 • 10~ 3 Z Q at z — 11 which is the redshift of 
reionization derived from the WMAP 3 year Compton optical depth (Spergel et al. 2006). 
This value is comfortably larger than the minimum metallicity required for reionization (by 
PopIII stars, Z « 1.2 • 10~ 4 Z Q , Stiavelli et al. 2004). This extrapolation could thus be 
considered as conservative. 

An additional test can be performed in terms of the predictions given by models for the 
formation and chemical enrichment of DLA systems based on cosmological hydro-dynamical 
simulations (Cen et al. 2003) or on semi-analytical prescriptions (Johansson & Efstathiou 
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2006). These models confirm a relatively slow but progressive decrease of the metallicity, 
while they evidence a sharp drop in the comoving density of neutral gas in DLA systems at 
z > 7 (Johansson & Efstathiou 2006, private communication), that however could be due to 
the built-in assumptions, e.g. coeval evolution with equal ages. 

To check what would be the consequences for our estimates in a scenario where the 
number of absorbers drops significantly at z > 7, we have run some Monte Carlo simulations 
assuming that there are no absorbers at z > 7 and compared the obscuration given by 
these models with that predicted by those with extrapolation at constant A. The results - 
presented in detail in Sec. 5 - evidence a modest variation in the average transmission. This 
is easily understood as, given the exponential decrease of k, the impact of z > 7 absorbers 
is limited. 



In order to compute the probability distribution for the optical depth we resort to a 
Monte Carlo code. Our code accepts general input functions <&(Nd), n(z), k(z) and £(A) and 
generates the chosen number of discrete realization for the DLA system distribution along 
lines of sight for a given emission redshift z e . The cumulative probability distribution for 
the optical depth to redshift z e is then built. 

For each discrete realization, we begin by integrating up to z e the redshift distribution 
of DLA systems n(z) so as to obtain the expected total number of absorbers along the line 
of sight: 



A Poisson random variable with mean ritot representing the realized number of absorbers 
is generated using a standard subroutine from Press et al. (1992). The redshift for each 
absorber is then randomly assigned by inversion of the primitive of n(z), evaluated numer- 
ically. Similarly the column density for each absorber is assigned via random sampling by 
inversion of the primitive for &(N d ). Once the redshift and column density for each absorber 
has been assigned, the optical depth at the observed wavelength of interest is computed for 
each absorber via Eq. 8. The total optical depth is given by summing over all the absorbers 
along the line of sight. 

The accuracy of the code is estimated by evaluating the variance for selected levels 
in the optical depth distribution (median, upper and lower 1,2a points). In addition the 
effective optical depth is compared to the analytical value from Eq. (9). The effective optical 
depth and the median are evaluated, at a fixed number of discrete realizations, with greater 
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Our Monte Carlo Code 




(16) 



- 11 - 



accuracy than the 1 and la points. Especially for the 2a contours, one order of magnitude 
more realizations are needed for an accuracy comparable to the one reached for the effective 
optical depth. For our purposes we are satisfied with an absolute error below 10~ 3 for r^ e ^\ 
This is reached with about 10 5 random lines of sight. Depending on the redshift of emission 
(lines of sight for low z e have an expected number of intervening DLA system much less than 
one, so the relative variance in the MC code is higher) this translates into a relative error 
below 1CT 2 at low z e and of about 2 — 3 • 1CT 3 at high z e for the MC simulations presented 
in this paper (see Fig. 2). 



5. Results: absorption probability distribution 

In Figs. 3-4 we plot the effective optical depth r^ 6 -^ (at the emitter rest-frame A e = 
0.14 \xm) to a given redshift, obtained using our models E101_a and E101_b, which employs 
the best fitting value for A from the CORALS survey (that represents a generous upper 
limit for A, see Sec. 3) and, respectively, a K = 1 (50% of metals in dust) for model "a" and 
a K = 0.5 (25% of metals in dust) for model "b". 

The effective optical depth peaks at about z ~ 5 with a maximum value below r^-^ < 
0.08 for the model with the highest dust-to-gas ratio. As the emission redshift increases 
the optical depth at fixed emitted wavelength decreases. In fact, despite the increase in the 
redshift density n(z), high z absorbers are characterized by a lower metallicity, decreasing 
exponentially with z in our model, while absorbers at lower z are traversed by light that 
has been redshifted toward progressively longer wavelengths, where the absorbers are more 
transparent. This explains the shape of the differential contribution to the effective optical 
depth shown in Fig. 5. Even for z > 10 observations the main contribution to absorption 
comes from systems at 2 < z < 5, which is a range probed with the highest accuracy by 
current DLA surveys. 

The optical depth distribution is characterized by a small number of highly obscured 
lines of sight, while the vast majority is almost dust-free: e.g., for our standard model 
(E101_a), along a random direction r < 0.1 with probability pa 0.8. Only 5% of the lines 
of sight may have r > 0.35 to z — 5, while to z = 20 we have r < 0.1 with probability 
0.95. In the E101_b model the optical depth is below 0.2 with probability 0.95 at z pa 5 and 
declines below 0.1 for z > 12 with probability 0.95. The optical depth distribution (Figs. 3-4) 
shows a sharp rise from to the value of the minimum optical depth for a single absorber at 
the redshift where the probability of having a clear line of sight falls below the probability 
associated to the line plotted. 
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Decreasing the dust-to-gas ratio (models E101_a-E101_c) leads to a corresponding quasi- 
linear decrease in the effective optical depth (see Eq. 9 and compare Fig. 4 - model E101_b 
- with Fig. 3 - model E101_a ). 

Even if we consider alternative models the expected effective optical depth does not 
change dramatically. In Fig. 6 we report the optical depth distribution for our models 
calibrated to the SDSS DLA survey data. The two different models considered (Pr05_P 
-power law- and Pr05_r -gamma function- for $(A^)), have negligible differences between 
each other in terms of the resulting effective optical depth. The SDSS data suggest an 
effective optical depth that is about 30% smaller than that from the CORALS data. 

In Fig. 7 we explore the effects of variations of different assumptions for the model- 
ing of the absorber distribution. One important parameter that is difficult to constraint 
observationally is the cut-off value N max for the power law form of ^A*^) used to fit the 
CORALS data. We have considered two additional models (E101_d and E101_e) with cutoffs 
two and ten times higher than E101_a (see Tab. 1). The effective optical depth, shown in 
Fig. 7, changes by At^-^ < 0.02 going from E101_a to E101_e. Absorbers with high dust 
column densities, that may have been missed in the CORALS survey due to small number 
statistics (and that are likely to be missed in optically selected surveys like SDSS), have only 
a modest effect on the expected average transmission. One caveat is that this conclusion 
has been reached by extrapolating a power law column density distribution into a region of 
the parameter space (Nj > 10 22 cm~ 2 ) where there are no observational constraints. These 
systems may well be a distinct population of absorbers with column density and dust-to-gas 
ratio distributions different from those of the observed DLAs. In particular an absorber with 
N d > 10 22 cm~ 2 may contain H 2 and hence have an higher dust-to-gas ratio. An upper limit 
to the uncertainty introduced on the expected average transmission by a hypothetic popula- 
tion of "bricks" absorbers can be estimated from the CORALS survey. At 95% of confidence 
level this population influences less than 2.9% of the lines of sight. At this confidence level 
the maximum displacement introduced in r^ 6 -^ is 

& T (eff) < o.03, in good agreement with 
the estimate Ar* 6 -^ < 0.02 that we obtain with our E101_e model (for the effects of dust-rich 
absorbers see also the model E101_h discussed below). 

Model E101_f, shown in Fig. 7, is characterized by a cut-off of the absorbers distribution 
at z — 7, while otherwise coincides with the model E101_a. This allows to quantify the 
uncertainties associated with the extrapolation of our models to redshifts where DLA data 
are unavailable. The possibility that the DLA number density may drop significantly at 
high redshift is hinted by semi-analytical models (Johansson & Efstathiou 2006). Even in 
the extreme scenario of model E101_f, where the DLA number is set to for z > 7, the 
difference in the average effective optical depth is only Ar^ < 0.008. 
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Model E101_g investigates the effects of the presence of a sub-population of DLA systems 
at z < 1.8 with high metallicity (Wild & Hewett 2005; Wild et al. 2006). This model has 
been constructed starting from the standard E101_a and assuming that one third of the DLA 
systems at z < 1.8 have solar metallicity. The presence of a significant number of these 
systems at higher z is unlikely, as the metallicity in the CORALS survey is significantly 
sub-solar and consistent with that measured in optically selected surveys (Akerman et al. 
2005). The effect of this population of absorbers is to enhance the optical depth up to z < 6 
(see Fig. 7). However as the emission redshift further increases their influence is significantly 
reduced and becomes negligible for z > 10. 

Model E101_h continues to investigate the effects of a small population of optically thick 
absorbers (like the inner regions of spiral galaxies or dusty star bust galaxies like M82) 
that could have been missed in the CORALS survey due to the limited number of lines of 
sight probed. We assume that a random line of sight intersects a number of galaxies drawn 
from a Poisson distribution with average 0.025 1 . We assume that each galaxy introduces an 
optical depth of 0.5 (estimated from Holwerda et al. 2005). The results of the Monte Carlo 
simulation, shown in Fig. 8, indicate that the effective optical depth is slightly higher in 
this case (i.e. the average transmission is marginally lower with respect to E101_a). At the 
level of the optical depth distribution, only the contour lines associated to transmission <C 1 
are influenced, i.e. those directly affected by lines of sight intersecting a galaxy. Differences 
between the E101_h and the E101_a models appear significant only for the top 10 % of the 
distribution. 



5.1. Reddening 

The shape of the effective extinction curve _E[£(A)] = log E[q(\, z e )]/ log E[q(\B, z e )] is 
only marginally dependent on the model used or on the emission redshift considered (see 
Fig. 9). We can empirically fit in the range A e G [O.l^m, 8/j,m] the effective average extinction 



1 This value for the covering factor has been estimated as sum of two contributions at high and low 
redshift. For the high redshift contribution we have analyzed the Hubble Ultra Deep Field. In the 2775 band 
the covering factor of pixels brighter than nii = 33 is w 0.01. This number is broadly compatible with what 
is derived from the luminosity function of Lyman Break galaxies in the redshift range 2 < z < 6 (Steidcl 
et al. 1999; Ferguson et al. 2004; Bouwens et al. 2005) and extrapolated down to z = 1. To estimate the 
covering factor due to galaxies at z < 1 we have considered four SDSS random fields in the i band for a total 
area of « 0.15 deg 2 . For each field we have applied a cut at +2.5ct from the sky level and then removed 
isolated pixels. This analysis provides an estimate of the covering factor of w 0.015, with single field values 
in the range [0.008; 0.022]. 
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curve for emitters at redshift z > 1 with a simple power law in the form: 

«A e )]=^-(£) C , (17) 

with i] = —0.55, k = 0.755 /im and ( = 0.87. 

With this "universal" extinction curve, we can infer the typical values for the average 
reddening, e.g. E(B — V), that can be estimated as a fixed fraction of the effective optical 
depth r* 6 ^ at a reference wavelength. 

6. Discussion 

We present a model for quantifying the effects of absorption due to dust in DLA systems 
along the line of sight for sources up to z = 20. 

The effective optical depth to a given redshift as a function of the emitted frequency A e 
can be evaluated analytically by using Eq. (9). This allows to obtain immediately an order 
of magnitude estimate of the effects of the dust absorption on the average transmission for 
the class of observations one is interested in. For a better characterization of the effects of 
the extinction, we study by means of a Monte Carlo method the distribution of the optical 
depths to a given redshift, setting upper and lower limits on the dust extinction. 

Under the reference scenario, that accounts for a large fraction of metals in dust grains 
(model E101_a with 50% of the metals in dust) the effects of dust obscuration remain modest 
even for very high redshift, with an optical depth at A e = 0.14/im below 0.4 with probability 
0.95 for z e pa 3. As the emission redshift increases the optical depth decreases and for z e > 15 
our modeling predict r < 0.1 with probability 0.95. We have explored several alternative 
possibilities for the input parameters finding that the effective optical depth varies within a 
factor 2 at most even when a population of optically thick absorbers like the central parts of 
star-forming galaxies is taken into account. Therefore the loss of sensitivity and the effects 
of reddening are not expected to significantly influence high-z observations with the James 
Webb Space Telescope. 

In the future we plan to extend the present framework to include additional effects on 
the transmission along the line of sight, such as gravitational lensing magnification and de- 
magnification, that may be significant for explaining the observed number counts of bright 
quasars in the SDSS Damped Lyman a survey (Prochaska et al. 2005). 

It is our pleasure to thank Mike Fall, Harry Ferguson, Benne Holwerda and Jason Xavier 
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for providing additional unpublished data from his semi-analytical model for formation and 
evolution of DLA systems. We are grateful to the referee for constructive suggestions that 
have improved the paper. This work was supported in part by NASA JWST IDS grant 
NAG5-12458. 
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Fig. 1. — Maximum likelihood estimation for the parameters A (left panel) and a2 (right 
panel), based on the data from Ellison et al. (2001), dotted line, and from Prochaska et al. 
(2005), solid line. For each dataset we plot the likelihood L (log scale) of the parameters 
normalized to the maximum value. 
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Fig. 2. — Absolute (left) and relative (right) error for the effective optical depth computed 
via our MC method using 4 • 10 5 random lines of sight and compared with the analytical 
value from Eq. 9 for the models E101_a, E101_b and E101_c. 
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Fig. 3. — Left Panel: Effective optical depth r^ 6 -^ as a function of redshift at A e = 0.14/im 
(bold line, solid) for the E101_a model. The green dashed-dotted lines represent, starting 
from the top in the left panel, the upper 95% contour in the distribution of optical depth, 
the upper 68%, the median, the lower 32% and the lower 5%. Right Panel: Cumulative 
probability distribution for the optical depth at A e = 0.14/xm along lines of sight to different 
redshifts (solid z = 5, dotted z = 10, dashed z = 15). The curves have been generated with 
a MC code using 4 • 10 5 random lines of sight. 
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Fig. 4. — Effective optical depth (left) and cumulative probability distribution for the optical 
depth (right) as in Fig. 3, for the E101_b model. 
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Fig. 5. — Contribution to the integral in Eq. (9) divided in redshift bins. The main con- 
tribution to the effective optical depth to high redshift sources is given by absorbers at 
2<z< 5. 
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Fig. 6. — Effective optical depth at an emitted wavelength A e = 0.14/xm for the two Pr05 
models compared with the E101_a model. 
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Fig. 7. — Effective opticai depth at an emitted wavelength A e = 0.14/im for different models. 
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Fig. 8. — Optical depth vs. redshift for the model E101_h compared with the model E101_a. 
Effective optical depth is the bold black line for E101_h and a bold dotted line for E101_a. 
The green dashed-dotted lines represent for E101_h, starting from the top, the upper 95% 
contour in the distribution of optical depth, the upper 68%, the median, the lower 32% and 
the lower 5%. The red dashed lines represent the same quantities for E101_a. 
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Fig. 9. — Effective mean extinction curve (£(A)) = ln(E[q(\)]/ln(E[q(\B)])) for the E101_a 
model at z = 2.5 (dashed) and z — 5, (dotted-dashed) compared with the input SMC 
extinction curve (solid). Our analytical fitting formula, Eq. 17 (dotted line), is shown for 
comparison. 
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Table 1: Adopted parameters 



ID 


A 


Nj 




1 ' mm 


1 "max 


a K 


Pr05T 


0.0715 


3.03 


1.8 


0.2 


OO 


1 


ID 


A 






N ■ 


N 


a K 


Pr05_P 


0.0715 




2.2 


0.2 


10 


1 


E101_a 


0.0910 




2.1 


0.2 


10 


1 


E101_b 


0.0910 




2.1 


0.2 


10 


0.5 


E101_c 


0.0910 




2.1 


0.2 


10 


0.75 


E101_d 


0.0910 




2.1 


0.2 


20 


1 


E101_e 


0.0910 




2.1 


0.2 


100 


1 


E101J 


0.0910 




2.1 


0.2 


10 


6(7 -z) 


E101_g 


0.0910 




2.1 


0.2 


10 


l+'DLAbricks' 


ElOlJi 


0.0910 




2.1 


0.2 


10 


l+'disks' 



Note. - 



Summary table with the parameters used to compute dust absorption. 



